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Summary 

We report 15 new avian influenza virus A/H5N1 haemagglutinin (HA) 
sequences sampled from visibly sick domestic poultry in southern Vietnam, 
between 1 January 2010 and 6 March 2010. These HA sequences form a new 
sub-clade of the clade 1 H5N1 viruses that have been circulating in Vietnam 
since 2003/2004. The viruses are characterized by a change from isoleucine to 
valine at position 514 (15 14V) and are 1.8% divergent at the nucleotide level 
from HA sequences sampled in Vietnam in 2007. Five new amino acid changes 
were observed at previously identified antigenic sites, and three were located 
within structural elements of the receptor-binding domain. One new mutation 
removed a potential N-linked glycosylation site, and a methionine insertion 
was observed in one virus at the polybasic cleavage site. Five of these viruses 
were sampled from farms where poultry were vaccinated against H5N1, but 
there was no association between observed amino acid changes and flock vacci- 
nation status. Despite the current lack of evidence for antigenic drift or 
immune escape in Vietnamese H5N1 viruses, continued surveillance remains a 
high priority. 



Over 500 human cases of H5N1 influenza virus infection 
have been recorded since the emergence of this highly 
pathogenic virus subtype in 1996/1997 (Peiris et al., 2007; 
Tarantola et al., 2010). Although the global incidence of 
human H5N1 cases has declined recently, surveillance and 
molecular characterization of highly pathogenic avian 
influenza (HPAI) viruses remain a high priority. The 
majority of human cases have occurred in three countries, 
Indonesia (174), Egypt (130) and Vietnam (119); all of 



which have high rates of household poultry ownership, a 
large network of markets with high levels of poultry sales 
and consumption (Domenech et al, 2009) and frequent 
outbreaks and high prevalence of HPAI in domestic poul- 
try. Vietnam is home to more than 220 million domestic 
poultry, with more than 40% of households owning 
chickens or ducks. H5N1 avian influenza has been a con- 
stant threat as a potentially fatal zoonotic infection, caus- 
ing regular outbreaks between December and February 
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every year. Although the majority of human H5N1 cases 
in Vietnam occurred in 2004 and 2005, the past 2 years 
have seen 12 cases including seven fatalities. 

The global phylogenetic structure of avian H5N1 
viruses is described by WHO -defined clades and sub- 
clades based on haemagglutinin (HA) sequences (WHO/ 
OIE/FAO H5N1 Evolution Working Group, 2009). Viet- 
namese H5N1 viruses comprise clade 1, clade 2.3.2 and 
clade 2.3.4 (Le et al., 2008; Nguyen et al., 2008). Clade 7 
viruses were introduced across the Vietnamese-Chinese 
border in 2008 (Nguyen et al., 2009), but have not yet 
become established in Vietnam. To date, H5N1 circula- 
tion in southern Vietnam has been primarily, but not 
exclusively, limited to clade 1 viruses (Wan et al., 2008). 
The genetic structure of H5N1 virus populations in Viet- 
nam suggests that gene flow is limited between northern 
and southern Vietnam (Smith et al., 2006a). 

Since October 2005, the Vietnamese government has 
pursued ambitious biannual mass vaccination campaigns 
for domestic poultry (Hinrichs et al., 2010). This vaccina- 
tion strategy is one of the important explanatory variables 
of the reduced frequency of outbreaks and lower overall 
levels of infection (Domenech et al., 2009; Henning et al., 
2009). Approximately 250-350 million poultry are vacci- 
nated in March/ April and August/September each year 
with an inactivated vaccine generated by reverse genetics; 
the current formulation (Re-1) contains the HA and 
neuraminidase of A/Goose/Guangdong/ 1/96. Re-1 is now 
only distantly related to currently circulating viruses; 
however, it provides good protection against contempo- 
rary clade 1, 2.2 and 2.3.4 viruses (Tian et al., 2010). In 
2011, Vietnam plans to update the vaccine strain to Re -5 
( A/duck/ Anhui/ 1/2006), a clade 2.3.4 virus more closely 
related to Chinese circulating strains (Chen, 2009). Moni- 
toring genetic diversity and ongoing evolution in vacci- 
nated populations is critical for identifying changes in 
virulence and immune escape, but surveillance is compro- 
mised by the fact that some chickens seroconvert without 
showing signs of disease (Henning et al., 2010). One 
explanatory hypothesis for the lack of clinical signs posits 
the circulation of H5N1 variants exhibiting low pathoge- 
nicity in chickens; however, virological evidence to con- 
firm this is not yet available. 

Methods 

We sequenced the HA of H5N1 viruses obtained from 
diseased domestic poultry during passive surveillance 
activities of January-March 2010 in four provinces of 
southern Vietnam: Khanh Hoa, Ben Tre, Soc Trang and 
Ca Mau. 

Sampling was conducted by veterinary technicians from 
the respective provincial sub -departments of animal 
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health as per national guidelines in response to suspected 
HPAI outbreaks. Cloacal specimens were collected from 
one to three birds per flock and screened by RT-PCR 
using OIE standard protocols (OIE, 2008). Sixteen farms 
reported diseased poultry during the 3 -month study 
period, and 17 samples collected from these farms tested 
H5 -positive. Four farms reported a history of vaccination 
within the previous 3 months. Direct sequencing of HA 
was attempted on all samples as previously described 
(Guan et al., 2002), and 15 HA sequences were success- 
fully generated (9 from meat ducks, 4 from chickens and 
one each from a quail and Muscovy duck; Genbank 
accession numbers CY081026-CY081040). 

To place the 15 new HA sequences in regional context, 
we downloaded all HA sequences (>1600 nt) of Vietnam- 
ese avian H5N1 viruses from GenBank (n - 105; 33 
chicken, 71 duck and 1 quail), only including sequences 
from the Influenza Genome Sequencing Project (IGSP) 
(Ghedin et al., 2005) as the quality control and sequenc- 
ing redundancy in the IGSP makes it more likely that 
these sequences are free of sample contamination and 
sequencing errors (Boni et al, 2010). Our 15 HA 
sequences were aligned with the 105 GenBank sequences, 
16 clade markers, three clade 0 sequences as outgroups 
and three vaccine strains used in China (Chen, 2009), 
which are the most likely vaccine candidates to be consid- 
ered for use in Vietnam. Alignments were performed with 
MUSCLE v3.8 (Edgar, 2004), and maximum-likelihood 
phylogenetic inference was performed with RAxML using 
1000 bootstrap replicates (Stamatakis, 2006; Stamatakis 
et al., 2008). Figure 1 shows the phylogeny of all 142 
sequences. 

To investigate whether sequences outside Vietnam may 
be related to the 2010 strains sequenced in this study, we 
used a broader set of sequences to infer a second phyloge- 
netic tree using both IGSP and non-IGSP sequences. We 
downloaded all subtype H5N1 HA sequences (>1600 nt) 
from Cambodia (n = 25), Laos (n = 36), Vietnam 
(n = 293), the two Chinese provinces of Yunnan (n - 68) 
and Guangxi (n = 78) that border Vietnam, and the Chi- 
nese province of Guangdong (n = 37) that is close to 
Vietnam and an important epicentre of avian H5N1 cir- 
culation. Chinese sequences sampled before 2003 were 
excluded to reduce the size of the data set. These 
sequences were aligned with clade markers, vaccine strains 
and outgroups as above. Figure 2 shows a maximum- 
likelihood phylogenetic tree (RAxML, 100 bootstrap 
replicates) of these 574 sequences. 

Results and Discussion 

The fifteen 2010 sequences formed a clear bootstrap-sup- 
ported sub-clade within the clade 1 sequences, with a 
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Fig. 1. Maximum-likelihood phylogenetic tree of 142 Vietnamese avian H5N1 HA sequences. Tree rooted on Gs/Gd/1/96. Only sequences avail- 
able through the Influenza Genome Sequencing Project (National Institutes of Health, USA) were included in the analysis. Bootstrap values >70% 
are shown on the branches. Sequences from 2010 are colour-coded by location: Khanh Hoa (red), Ben Tre (orange), Soc Trang (green) and Ca 
Mau (blue). Grey boxes on branches denote the five sequences that were sampled from vaccinated farms. Vaccine strains used in China (Chen, 
2009) shown in magenta. Model of evolution used by RAxML: GTR-r, a = 0.422679, AC = 1.149255, AG = 5.429345, AT = 0.628623, 
CG = 0.232736, CT = 7.840430, GT = 1 .0. 
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Table 1. 15 HA sequences sampled during January-March 2010 in southern Vietnam 
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Sequence name 


Date of collection 


Location (province) 


New amino acid changes from clade 1 2007 sequences 


A/ChickenA/ietnam/1/201 0* 


Jan 1 


Ca Mau 


1514V, A214S 


A/DuckA/ietnam/1/201 0 


Jan 7 


Ca Mau 


1514V, A185E, N309T, D45G, R458K 


A/DuckA/ietnam/2/201 0* 


Jan 10 


Ca Mau 


1514V 


A/DuckA/ietnam/3/201 0 


Jan 17 


Ca Mau 


1514V, S155G, N236S; insertion (-)326M 


A/DuckA/ietnam/4/201 0 f 


Jan 24 


Ca Mau 


T204M 


A/DuckA/ietnam/5/201 0 f 


Jan 25 


Soc Trang 


S155N 


A/ChickenA/ietnam/2/201 0 


Feb 2 


Soc Trang 


1514V, N309T 


A/DuckA/ietnam/8/201 0 


Feb 12 


Soc Trang 


1514V, A185E, V200I 


A/DuckA/ietnam/7/2010 


Feb 17 


Soc Trang 


1514V, N309T, A184G 


A/ChickenA/ietnam/4/201 0 


Feb 23 


Khanh Hoa 


1514V 


A/QuailA/ietnam/1/2010 


Feb 23 


Khanh Hoa 


1514V, C42W, K161R 


A/Chicken/Vietnam/3/201 0 


Feb 25 


Khanh Hoa 


1514V 


A/Chicken/Vietnam/5/201 0* 


Feb 25 


Khanh Hoa 


1514V, 171 T 


A/DuckA/ietnam/9/201 0** 


Mar 6 


Ben Tre 


1514V, T36A, P136S, L41H, C42R 


A/DuckA/ietnam/1 0/201 0* :i: 


Mar 6 


Ben Tre 


1514V, T36A, P136S 



*vaccinated farm. 

t Not sequenced through position 514. 
Samples from the same farm. 

Boldface amino acid changes observed in more than one virus. 



mean distance of 31 nt (1.8%) from the ancestral 2007 
viruses (Figure 1). The 2010 clade is characterized by a 
substitution from isoleucine to valine at position 514 
(15 14V; H5 numbering as in WHO Global Influenza Pro- 
gram (2005) and Smith et al. (2006a)), with four new 
amino acid changes observed in at least two sequences 
(Table 1). Two previously polymorphic sites (K48R and 
N168D) exhibited dramatic increases in frequency in the 
2010 clade as compared to the 2007 clade. Five new 
amino acid changes were located in previously identified 
antigenic sites (T36A, D45G, P136S, S155N/G and 
A185E) (WHO Global Influenza Program, 2005; Stevens 
et al, 2006). Three amino acid changes (A184G, A185E 
and A214S) occurred within predicted loop structures of 
the receptor-binding domain (Stevens et al, 2006), and 
the S155N/G mutations remove a potential N-linked gly- 
cosylation site (Smith et al, 2006a; Yen et al., 2009). One 
virus had a single methionine insertion at the polybasic 
cleavage site. The phenotypic effects of these changes will 
be characterized in future analyses. 

When considered in a broader phylogenetic context - 
including IGSP and non-IGSP sequences from Vietnam, 
Cambodia, Laos and southern China - the 2010 Vietnam- 
ese HA sequences still form a bootstrap-supported (99%) 
sub-clade of the clade 1 viruses. Although it is likely that 
clade 1 viruses circulating in 2008 and 2009 were geneti- 
cally close to the 2010 viruses described here, no Viet- 
namese H5N1 sequences from 2008 or 2009 are available 
in GenBank. 

The 2010 viruses sampled in this report originated 
from chicken and duck farms with flock sizes ranging 



from 4 to >8000 birds. The outbreak farms experienced 
significant losses, with most farms experiencing between 
10% and 45% flock mortality. Four of the 15 farms 
reported flock vaccination within the previous 3 months; 
however, we do not have access to more precise informa- 
tion or flock histories to document what percentage of 
the flock was vaccinated nor to confirm whether all ducks 
received the recommended two doses. It should be noted 
that even if sera had been collected from the outbreak 
farms, it would not have been possible to confirm vacci- 
nation status, because the current vaccine formulation in 
Vietnam (inactivated Re-1 unmarked virus) generates an 
immune response that is indistinguishable from natural 
infection with field virus. We strongly recommend that 
future outbreak response activities more precisely docu- 
ment the details of flock vaccination. 

Many factors on both an individual and flock level may 
affect vaccine efficacy, such as problems with cold chain 
or poor vaccine quality, improper dosing (inappropriate 
age, failure to administer booster) or population turnover 
within the flock. Previous studies indicate that vaccination 
under field conditions may fail to achieve the same magni- 
tude of immunity as observed in experimental validations 
(Henning et al., 2010). In particular, seroconversion is 
known to be particularly poor in ducks and requires a 
minimum of two doses, recommended at 7 and 21 days, 
for full protection (Henning et al., 2010). Thus, it is diffi- 
cult to determine the true immune status and susceptibil- 
ity to infection or outbreaks of any particular flock. 

Our results concur with other reports suggesting that 
current vaccination campaigns are failing to provide 
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adequate protection from infection and disease in vacci- 
nated flocks and may thus be generating partially immune 
populations that could select for new antigenic variants 
of H5N1 (Smith et al., 2006b; Domenech et al, 2009; 
Henning et al., 2010). However, there is as yet no evidence 
that vaccination failure is the result of viral immune 
escape. Complete antigenic characterizations of the viruses 
sequenced in this study will be required to fully evaluate 
the phenotypic implications of the observed mutations. In 
addition, it must be emphasized that the presence of new 
amino acid changes described in the 2010 H5N1 Vietnam- 
ese viruses did not correlate with vaccination status; two 
of the viruses from vaccinated farms did not show new 
amino acid changes at antigenic sites. 

The combination of aggressive vaccination campaigns 
and HPAI outbreak control measures in Vietnam has 
substantially reduced H5N1 virus circulation throughout 
Vietnam since the first outbreaks in 2003-2005. However, 
HPAI continues to circulate among poorly vaccinated 
flocks. As H5 avian viruses have demonstrated capacity 
for immune-escape evolution in the past (Lee et al., 2004; 
Domenech et al., 2009), it is clear that sustained and 
vigilant surveillance of field viruses remains critically 
important. 
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